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Dynamic spin correlation function of the S = 1/2 antiferromagnetic spin chain with alternating 
nearest neighbor and competing next nearest neighbor interactions is studied using the numer- 
ical diagonalization method. We find characteristics of the correlation function, which are the 
consequence of the dimer order in the ground state: (1) A typical singlet-triplet gap is observed 
and such a lowest-energy triplet state is found to form a bound state, showing strong peaks 
below a continuum band with weak but distinct upper and lower boundaries. (2) The intensity 
of the bound state shows a quite different wavevector-dependence from that of the corresponding 
intensity at the lower boundary in the usual antiferromagnet. These features are examined by 
the variational calculation and are favorably discussed in connection with recent experimental 
results of neutron inelastic scattering in the inorganic spin-Peierls system, CuGeOs. 
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§1. Introduction 

It is realized recently that a spin gap due to the sin- 
glet spin dimer is rather common features in low dimen- 
sional quantum spin systems and plays important roles in 
the mechanism of the oxide superconductor. The spin- 
Peierls order, which is a consequence of the combined 
effects of the lattice distortion and the singlet dimer, is 
one of the highlights in this field and has attracted a 
renewed interest since the discovery of the inorganic Cu 
compound, CuGe03.EP 

Because of large and high-quality samples, it becomes 
able to perform neutron inelastic scattering (NIS) for 
CuGeC>3, 2,3 ' which yields invaluable information on the 
dynamic properties of the spin-Peierls system. We be- 
lieve that this is the most important progress in exper- 
iments comparing with .the situation in organic materi- 
als about 15 years agoo^ According to the experimen- 
tal results, the dynamic spin correlation function seems 
to be, at least at a first glance, very similar to that of 
the typical one-dimensional antiferromagnet. It is, how- 
ever, realized several important differences: (1) At the 
wavevectors Q — and 7r in units oil/ a where a is the 
lattice constant in the uniform phase, it exhibits a typi- 
cal gap. (2) In the whole wavevector space, bound state 
peaks seem to exist below the continuum band which 
has a rather distinct lower boundary.© (3) The intensity 
of the bound state decreases with decreasing Q from 7r 
but the decreasing rate is much slower than that in the 
typical one-dimensional antiferromagnet. a 

It is the purpose of this paper to reveal characteris- 
tics of the dynamic spin correlation function of the one- 
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dimensional quantum spin chain which has the dimer or- 
der in the ground state.© We do not aim, at present, at 
making a quantitative comparison with the experimental 
results but want to clarify the characteristic features of 
the correlation function of the system described by the 
following Hamiltonian: 

n = O 1 + (- 1 )'^) S « ■ + jS„ ■ S n+2 ], (1.1) 

n 

where S„ represents the spin-1/2 operator at the nth 
site. 5(1 > 5 > 0) is the parameter describing the 
bond-alternation and j(> 0) is the strength of the next 
nearest neighbor (nnn) interaction divided by that of 
the nearest neighbor (nn) interaction which is set to be 
unity. This Hamiltonian is believed to be appropriate 
for CuGeOs, 6 ' 7 ' although the parameter values for it are 
still controversial. 

The ground-state phase diagram of the system de- 
scribed by the Hamiltonian (1.1) is shown in Fig. I, 8 ' 9 ) 
where the space is spanned by the two parameters, 5 
and j, both of which lead cooperatively to the dimer 
ground state with an excitation gap. As is seen in Fig. 
1, the dimer order (D) takes place everywhere except 
for the 6 = and j < j c (— 0.24) 10,11 ) line, where the 
ground state is the usual gapless spin fluid (SF) state 
as in the case of S = and j = 0. On the disorder 
line, 5 = — 2j + 1, which is shown by the dashed line 
in Fig. 1, the ground state is the direct product of the 
singlet dimers of nn spins so that a spin correlates only 
with the other spin in the dimer. As the ground state 
on this disorder line is a prototype of the dimer order, 
we focus our attention in this paper mainly on the dy- 
namic spin correlation function in the case of the dimer 
ground state on this disorder line, which starts from the 
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independent dimer point, (<5 = l,j = 0), and ends at 
the point, (S — 0,j = 0.5), called the Majumdar-Ghosh 
(MG) point. It is worth noticing that the dimer ground 
states are doubly degenerated on the line, 6 = and 
j > j c , due to the translational symmetry of the Hamil- 
tonian while a finite value of 5 lifts the degeneracy. This 
observation leads us to the fact that there is a domain 
wall excitation in the former case but not in the latter 
case, as will be discussed in the following sections. 

In the next section, we describe our numerical method 
to reproduce the dynamic spin correlation function and 
present the result with some discussions. The variational 
calculation is introduced in §3 to extract a physical pic- 
ture from the numerical result. The last section is de- 
voted to our conclusion. 

§2. Numerical Calculation 

In this section, we first summarize briefly the—numeri- 
cal method developed by Gagliano and Balseirojij) which 
enables us to calculate the dynamic spin correlation func- 
tion S zz (Q,lo), making this paper self-contained. Then, 
the result of S zz (Q,ui) will be presented in a graphical 
form and its characteristics will be discussed. 

We define S zz (Q,uj) by 

i r°° 

S zz (Q,u;) = — / die— 4 <0|^ Q S^)|0)(2.1) 
= — lim 

x hn[(0\Si Q {u -ie + E - W) _1 S$|0>], (2.2) 
where 

S Z Q = N- 1 / 2 Y / cMiQn)Sn, (2-3) 

n 

S^(t) = exp(mt)S^exp(-mt), (2.4) 

Im[- ■ •] means the imaginary part and |0) denotes the 
ground state, which is obtained using the numerical di- 
agonalization method with the standard Lanczos algo- 
rithm. Once the ground state of a finite system is ob- 
tained, then S zz (Q,uj) can be calculated using the con- 
tinued fraction expansion (CFE) for the resolvent in the 
tridiagonal basis of the Hamiltonian with the initial vec- 
tor Sq\Q). The tridiagonal expansion of the Hamiltonian 
can be carried out numerically again with the Lanczos 
algorithm. Then, we obtain spectral weighs and poles 
of the S zz (Q,u>) with the CFE representation up to the 
order of 70. The convergence of the spectral weights 
and the pole positions are rather rapid in the low energy 
excitations, while it is worse in high energy excitations 
although they have only weak spectral weights. 

In Fig. 2, we present our numerical results of the dy- 
namic correlation function S zz (Q,u>) for N = 24 for sev- 
eral parameter values on the disorder line together with 
the result of an exceptional case, the antiferromagnetic 
chain (j = and (5 = 0), for comparison. 

In Fig. 2 (a), we show a typical S zz (Q,ui) for the 
nn Heisenberg antiferromagnet. The spectrum consists 
of the continuum band with the lower and the upper 
boundaries, which coincide with the known curves, lul = 



(7r/2)sin(Q) and wu = 7rsin(Q/2), respectively. We note 
that the intensity is distributed within the band although 
the intensity at the lower boundary is high and seems to 
diverge especially at Q = tt. The result agrees with the 
one known for this system, showing the validity of our 
numerical calculation, although our spectra consist of the 
discrete poles due to the finiteness of the system size. 

Now we proceed in cases where the ground state has 
the dimer order. We especially focus our attention on 
the perfect dimer cases, which occur on the disorder line 
as was mentioned in §1. The first example is the inde- 
pendent dimer point at which S zz (Q,oj) shows the delta 
function shape at to = 2 without dispersion. It is easy to 
understand this since each isolated dimer has the singlet- 
triplet gap of 2 in our units. A bulk triplet state is ob- 
viously the state where one of N/2 nn singlet dimers in 
the ground state is excited to the triplet state. As 5 de- 
creases and j increases, the degeneracy of these states is 
lifted and these states result in forming the triplet band, 
since the dimer propergates and spreads out due to the 
weak nn and nnn interactions. We show in Fig. 2 (b) an 
example of S ZZ (Q, oj) for this category, where the excited 
triplet states have a little dispersion. It is interesting to 
note that the state which consists of a linear combina- 
tion of such states with the wavevector tt/2 is the exact 
eigenstate of the Hamittcinian as far as the parameters 
are on the disorder lineJi-f 

At another end of the disorder line, the MG point, 
it is known that the triplet states form a band whose 
elementary excitation is a_domain wall connecting two 
degenerate ground statesO We note that the situation 
occurs only if 6 = 0. In Fig. 2 (d) the bound state 
predicted can be seen in the vicinity of Q = tt/2 and 
seems to merge into the lower boundary of the continuum 
band at a certain wavevector. This observation is again 
consistent with the known result l3 

The most typical spectrum of the dynamical correla- 
tion function on the disorder line can be seen in Fig. 2 
(c) and it also seems to be appropriate for the discussion 
of NIS in CuGeOa. First of all, we want to note that 
the lowest pole for each wavevector is isolated from oth- 
ers. Although the situation is smeared by the finite-size 
effect, we confirm this statement by the analysis of the 
size dependence of the positions and the spectral weights 
of such poles. This suggests that the lowest poles form 
the bound state splitting off from the triplet continuum 
band, which has rather distinct lower and upper bound- 
aries. 

Last point we would like to mention in this section 
is the wavevector dependence of the intensities at the 
lowest energy poles wb, which is shown in Fig. 3. In this 
figure S zz (Q,u)b)/S zz (k,ub) for j = 0.3 is plotted as a 
function of Q. We note that the denominator S zz (ir, ojb) 
itself decreases rather rapidly when S is introduced and 
the intensity is distributed to other wavevector modes. 
In addition, the higher energy tail shrinks as 5 increases 
so that the intensity at the lowest energy mode seems to 
have a relatively high intensity. Especially, at Q = tt/2 
the width becomes zero, so that the intensity can be 
comparable with that at Q = tt. This behavior of the 
intensities at the lowest energy poles is consistent with 
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the experimental observation. 

In the next section, we develop the variational calcu- 
lation, emphasizing an important role of the local triplet 
dimer in the singlet dimer background, to make clear the 
nature of the triplet excitations studied in this section. 

§3. Variational Calculation 

In this section, we pursue a variational calculation for 
the excited triplet states, assuming the complete dimer 
singlet ground state, in order to gain insights into the 
nature of the triplet bound state as well as of the triplet 
continuum band, which have been shown in the previous 
section. 

Before going into the detailed calculation, we consider 
how to construct appropriate variational states. The 
most plausible explanation for the low energy excited 
states in the system seems to be the one based on the 
domain wall scenario like in the case of the MG point. EJ 
However, we note, when we introduce the alternation, 
the domain wall is no longer an elementary excitation of 
the system since the ground-state degeneracy is lifted by 
the alternation. We suppose the following two types of 
the states as our variational states: (1) For the triplet 
bound state, we adopt the state in which one of N/2 
singlet dimers in the dimer ground state is excited to 
the triplet dimer. (2) For the triplet band states, we 
consider the states in which two triplet dimers among 
N/2 — 2 singlet dimers couple into triplet. Note that the 
former state is nothing but an extreme state of the two 
domain wall states, where two domain walls are located 
on nn sites. It is clear that these states are the low- lying 
eigenstates in the independent dimer limit (8 = l,j = 0). 

Based on the observation above, we build up the vari- 
ational states. For convenience, we employ the two-spin 
states on a stronger bond of a nn pair as local bases of 
the Hamiltonian. Then, the following symbols are intro- 
duced for the singlet state and the three components of 
the triplet state: 

\'n) = 2 1 ^ 2 | Traln+1 — IraTn+l) ! 
I In) = I tntri+l) , 
|o„) = 2~ 1 / 2 \ Uin+l + Intn+l) , 
I — In) = I Inln+l) , 

where n denotes the site and is assumed to be the even 
number. For example, the dimer ground state is ex- 
pressed by |D) = | •••••••• •}. We choose the following 

orthonormal states as the variarional bases for the low 
energy triplet states with the wavevector Q: 

\a Q )^(2/N)^J2 eiQm \ a ^ > 

n 

K, q2 ) = C{ qi ,q 2 )YM {qim+q2n) -e^ n+ ^\\b m ^.2) 

m,n 

where Q = qi + q2, and C{q\ 1 q2) denotes the normaliza- 
tion constant, and 

\a m ) = | ••• • »l m •••■■), (3.3) 

\b m ,n) = !•••• »l m • • • • • o„ • • • • •) . (3.4) 



These \a m ) and |6 m ,n) are depicted in Fig. 4. It is no- 
table that two triplet states, |6 m>n ) and \b n _ rn ), are super- 
imposed antisymmetrically so as to produce the S = 1 
continuum states. In addition, we deal only with the 
triplet states with S* = ^2 n S* = 1, for the computa- 
tional convenience. 

Now, we deduce the effective Hamiltonian within the 
above restricted Hilbert space. To this end, we fur- 
ther take into account the effect of states consisting 
of three triplet dimers within the second order pertur- 
bation theory, as the number of the triplet dimers is 
not conserved in this system. In other words, we in- 
tegrate the process into the effective Hamiltonian so 
that the scattering matrix elements between the two 
triplet states, (b qi>q2 \H\b q ' q ' ), and the diagonal ele- 
ments, {b qitq2 \T~t\b qi , q2 ), are renormalized. Then, wediag- 
onalize numerically the effective Hamiltonian matrix and 
obtain the excitation spectrum and the dynamic struc- 
ture factor S" 1 (Q,u>) defined by 

1 f°° 

S+-(Q,u)^— j^dte- lut (0\SZ Q S+(t)\0), (3.5) 
where 

S±=S x Q ±iS y Q . (3.6) 

Here we note that (Q, to) = 2S!?(Q,uj) since the 
system has the rotational invariance.EB' 

It is worth noticing that, in our variational calcu- 
lation, the off-diagonal elements of the Hamiltonian 
(&gi,«2 l^l a q)' s vanish at Q = tt/2 so that |oq) becomes 
an exact eigenstate in this case. Reminding the fact 
that we assume the complete singlet dimer ground state, 
which is true if the parameters are on the disorder line, 
we see that our variational calculation yields an exact 
solution at Q = tt/2 and on the disorder line. It is also 
remarkable to note that eq. (3.5) suggests the spectral 
weights of 5" 1 (Q, u) depend on how much the state |<zq) 
is mixed in the excited states. In fact, a contribution of 
the state |oq) seems to dominate 5 H (Q,oj) as far as the 
ground state has the dimer order. 

Next, we discuss the excitation spectra which are 
shown by the three dashed curves in Fig. 2 (b) and 
(c). These are, hereafter, denoted as u>b(Q), ^l{Q), 
and u>\j(Q) in the order of their energy, i.e., u>b(Q) < 
wl(Q) < wtj(Q)- The excitation spectrum forms the con- 
tinuum band in the region u>l(Q) < lu < u>u(Q), while 
the lowest energy state is a bound state, whose dispersion 
curve u>b(Q) agrees well with the lowest energy poles of 
S zz {Q,lj) for 5 > 0.2. 

It is quite reasonable that, the dispersion of all modes 
continuously tends to vanish with approaching the in- 
dependent dimer point 5=1 and j = 0, where the 
ground state consists of the N/2 independent dimers. 
In this limit, the first and the second excited states in- 
clude, respectively, one local triplet and two local triplets 
in the singlet dimer background. Considering these, we 
conclude that the continuum band consists of two local 
triplets: the upper and lower bounds are states where 
two triplets are located in the nn sites and in the distant 
sites, respectively. As is seen in Fig. 2 the dispersion 
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curve for the former coincides well with the numerical re- 
sult while that of the latter does not. We feel that, in or- 
der to obtain a reasonable result for the lower bound, the 
restriction of our variational calculation that the triplet 
is on nn sites must be released. In this sense our result 
becomes worse when S appraoches zero. 

As was mentioned above, in the limit of the MG point, 
(S = 0, j = 1/2), our perturbative approach fails to re- 
produce the continuum band obtained by the numerical 
diagonalization. For this case, the system possesses the 
full translational symmetry, giving rise to the degener- 
acy of the ground states. In other words, an elementary 
excitation in this case is not such a local triplet state 
but a spin-1/2 domain wall propagating on the lattice, 
as was shown by Shastry and Sutherland. EJ In fact, we 
show in Fig. 2 (d) their result for the upper and the 
lower bound of the continuum band with the dash-dotted 
curve. In this limit, we see the bound state in the vicin- 
ity of Q — 7r/2, which is the same as the one discussed 
above, although it merges into the continuum band at 
the wavevectors Q = 0.36-7T and 0.647T. 

In our variational approach, one additional bound 
state is observed in the energy spectrum, it lies between 
the lowest bound state and the lower boundary of the 
two triplet continuum as shown in Fig. 5. No corre- 
spondence is found in the result of the numerical diago- 
nalization. But it has only the quite small contribution 
to the scattering intensity. Perhaps it is caused by the 
strong restriction for the Hilbert space in our variational 
calculation, i.e. , each pair of spins that forms the lo- 
cal triplet state are bound more gradually in space, not 
bound strongly on the nearest neighbour sites. 

In conclusion, we see that our variational calculation 
gives us a reasonable explanation for the dispersion for 
the bound state as well as the continuum band obtained 
in the numerical calculation except for its lower bound- 
ary. Now, we turn to the wavevector dependence of their 
intensities on the lowest energy states. The behavior of 
the intensities coincides with that of the numerical re- 
sult, as shown in Fig. 6, not only qualitatively but also 
quantitatively, although this coincidence becomes worse 
when the parameters tend to the MG point. Decreasing 
S along the disorder line, we find the spectral weight of 
the lower energy mode decreases while that of the higher 
energy mode gets the intensity, except at Q = tt/2 where 
only the lower energy mode has the intensity. This ten- 
dency is again consistent with the result of the numerical 
calculation. 

§4. Conclusion 

In conclusion, we have found that in the S = 1/2 
antiferromagnetic chain with alternating nn interactions 
and competing nnn interactions, the dynamic correlation 
function shows characteristic features due to the dimcr 
order in the ground state: (1) The ground state is the 
disordered state so that the static structure factor, which 
is obtained by an integration of S zz (Q,lu) with respect 
to lo, has a relatively high intensity even at wavevec- 
tors other than Q = n. This is marked contrast to the 
case in the spin-fluid state. (2) The width of S zz (Q,ui) 
at Q = 7r/2 becomes narrow as the disorder line is ap- 



proached and eventually reaches zero at the disorder line. 
(3) The lowest pole for each wavevector splits off from 
the continuum band. We elucidated the nature of these 
states by the variational calculation assuming the triplet 
states based on the physical picture, and confirmed the 
facts mentioned above. 

All the behavior observed in our results is consistent 
with the experimental observation of NIS. 2,3 -' This sug- 
gests that the parameters of CuGe03 are close to those 
on the disorder line. We believe that the reason why 
the NIS result shows very different characteristic behav- 
ior from what we expect for a small dimerization in the 
nn spin-Peierls system is a manifestation of the impor- 
tant role played by the rather large competing nnn in- 
teraction. We have not aimed at comparing quantita- 
tively our theoretical result with the experimental re- 
sult of CuGe03 in this paper, because we feel that the 
following issues should be considered for a quantitative 
comparison. 

First we should throw light on the effect of the inter- 
chain interactions, since they are not so weak and said to 
be only one order of magnitude smaller than the intra- 
chain interactions. We believe, however, that the effect 
is reduced by singlet fluctuations caused by the rather 
large nnn interactions. Second a possibility that the 
dimerized lattice distortion can relax according to the 
spin state should be considered. If the lattice distortion 
follows completely the spin behavior, our Hamiltonian 
is inadequate in the sense that the distortion is treated 
as rigid. In such a case a soliton state like the one dis- 
cussed for the magnetic ground states in high fields 16 ' 17 -* 
may be a proper excited state. In this situation, infor- 
mation for the phonon is highly desired. In any case this 
is an interesting as well as a challenging problem for our 
future study. Lastly, the condition for the bound state 
to exist should be clarified from both experimental and 
theoretical sides. 

At the final stage of ouiistudy, we have been notified 
that Yokoyama and SaigaES have studied the same prob- 
lem using the same method. They claimed for CuGe03 
rather large value of j keeping a small value of 5, which 
is consistent with our conclusion in the sense that, in or- 
der to see characteristic features of the spin correlation 
function, the system should be close to the disorder line. 
We hope that these calculations stimulate further exper- 
imental work of NIS to elucidate the nature of the bound 
state as well as the continuum band, which enable us a 
quantitative comparison with theoretical calculations. 
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igure Captions 

1 Ground-state phase diagram in the parameter space 
spanned by the alternation parameter 5 and the 
competing next nearest neighbor interaction j. The 
solid line shows the spin-fluid (SF) phase without a 
spin gap. Other region is the dimer (D) phase with 
the spin gap. The disorder line is denoted by the 
dashed line and the Majumdar-Ghosh (MG) point 
is indicated. 



2 Dynamic spin correlation function of the finite spin 
system (N = 24), S zz (Q,u>), for representative val- 
ues of parameters, S and j. In (a) the dash-dotted 
curves show the rigorous upper and lower bound- 
aries of the triplet band, while in (b), (c) the approx- 
imate dispersion curves of the traveling triplet dimer 
in the singlet-dimer background are represented by 
the dashed curves. (See text.) In (d), the dash- 
dotted curves show the two-particle_scattering con- 
tinuum of the S = /2 domain wallal3' 



3 Intensity ratio of the lowest energy poles, 

S zz (Q , ujb) / S zz (n 7 uj b ), for j = 0.3 and representa- 
tive values of S as a function of the wavevector Q. 



4 Two types of triplet states adopted to construct the 
trial wave functions. In each state, the two arrows 
form the triplet state with S z + S z +1 = 1 while the 
solid rectangle denotes the triplet pair with S z + 
S z +1 = 0. The solid ellipse shows the singlet pair 
which is the same as in the ground state. 



Fig. 5 Triplet excitation spectra for j = 0.3 and S = 0.4. 
The solid lines represent the result obtained by the 
variational calculation while the open diamonds de- 
note the result obtained by the numerical diagonal- 
ization for the finite system (N — 24). The shaded 
area denotes the continuum band obtained by the 
variational calculation. 



Fig. 6 Wavevector dependence of the spectral weight on 
the lowest excited mode wb(Q)- The numerical re- 
sults of the finite chain (N — 24) for j — 0.1 and 
S = 0.8, j = 0.3 and S = 0.4, and j = 0.4 and 
S = 0.2 are depicted with the circles, the diamonds, 
and the triangles, respectively. The broken lines are 
to guide the reader's eye. The solid lines denote the 
variational results. 
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